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In this talk we present a novel framework that unifies the stunning success of MOND on 
galactic scales with the triumph of the ACDM model on cosmological scales. This is achieved 
through the rich and well-studied physics of superfluidity. The dark matter and MOND 
components have a common origin, representing different phases of a single underlying sub¬ 
stance. In galaxies, dark matter thermalizes and condenses to form a superfluid phase. The 
superfluid phonons couple to baryonic matter particles and mediate a MOND-like force. Our 
framework naturally distinguishes between galaxies (where MOND is successful) and galaxy 
clusters (where MOND is not): dark matter has a higher temperature in clusters, and hence 
is in a mixture of superfluid and normal phase. The rich and well-studied physics of su¬ 
perfluidity leads to a number of striking observational signatures, which we briefly discuss. 
Remarkably the critical temperature and equation of state of the dark matter superfluid are 
similar to those of known cold atom systems. Identifying a precise cold atom analogue would 
give important insights on the microphysical interactions underlying DM superfluidity. Tan- 
talizingly, it might open the possibility of simulating the properties and dynamics of galaxies 
in laboratory experiments. 


1 Introduction 

In the A-Cold-Dark-Matter (ACDM) standard model of cosmology, dark matter (DM) consists of 
collisionless particles. This model does exquisitely well at fitting a number of large-scale observa¬ 
tions, from the background expansion history to the cosmic microwave background anisotropies 
to the linear growth of cosmic structures 1 . 

On the scales of galaxies, however, the situation is murkier. A number of challenges have 
emerged for the standard ACDM model in recent years, as observations and numerical simula¬ 
tions of galaxies have improved in tandem. For starters, galaxies in our universe are surprisingly 
regular, exhibiting striking correlations among their physical properties. For instance, disc 
galaxies display a remarkably tight correlation between the total baryonic mass (stellar + gas) 
and the asymptotic rotational velocity, Mb ~ v^. This scaling relation, known as the Baryonic 
Tully-Fisher Relation (BTFR) 2,3 , is unexplained in the standard model. In order to reproduce 
the BTFR on average, simulations must finely adjust many parameters that model complex 
baryonic processes. Given the stochastic nature of these processes, the predicted scatter around 





Figure 1 Observed rotation curve for NGC1560 (blue points) 14 . The MOND curve (green) 15 offers a much 
better fit to the data than the ACDM curve (blue) 16 . Reproduced from 1 '. 

the BTFR is much larger than the observed tight correlation * * * 4 . 

Another suite of puzzles comes from the distribution of dwarf satellite galaxies around the 
Milky Way (MW) and Andromeda galaxies. The ACDM model predicts hundreds of small DM 
halos orbiting MW-like galaxies, which are in principle good homes for dwarf galaxies, yet only 
~ 20 — 30 dwarfs are observed around the MW and Andromeda. Recent attempts at matching 
the populations of simulated subhaloes and observed MW dwarf galaxies have revealed a “too 
big to fail” problem 5 : the most massive dark halos seen in the simulations are too dense to 
host the brightest MW satellites. Even more puzzling is the fact that the majority of the MW 6 
and Andromeda 7,8 satellites lie within vast planar structures and are co-rotating within these 
planes. (Phase-space correlated dwarfs have also been found around galaxies beyond the Local 
Group 9 .) This suggests that dwarf satellites did not form independently, as predicted by the 
standard model, but may have been created through an entirely different mechanism 6,10 . 

A radical alternative is Modified Newtonian Dynamics (MOND) 11,12 . MOND replaces DM 
with a modification to Newton’s gravitational force law that kicks in whenever the acceleration 
drops below a critical value ao- For large acceleration, a>ao, the force law recovers Newtonian 
gravity: a — on- At low acceleration, a <C a o, the force law is modified: a — y/a^ao. This simple 
empirical law has been remarkably successful at explaining a wide range of galactic phenomena 13 . 
In particular, asymptotically flat rotation curves and the BTFR are exact consequences of the 
force law. a MOND does exquisitely well at fitting detailed galactic rotation curves, as shown in 
Fig. 1. There is a single parameter, the critical acceleration ao, whose best-fit value is intriguingly 
of order the speed of light c times the Hubble constant Hq\ ao — \cHq ~ 1.2 x 1CT 8 cm/s 2 . 

However, the empirical success of MOND is limited to galaxies. The predicted X-ray temper¬ 
ature profile in massive clusters of galaxies is far from the observed approximately isothermal 
profile 18 . Relativistic extensions of MOND, e.g. 19 , fail to reproduce CMB anisotropies and 
large-scale clustering of galaxies 20 . The “Bullet” Cluster 21,22 , the aftermath of two colliding 
galaxy clusters, is also problematic for MOND 23 . 

2 Dark Matter Condensate 

In this talk, based on two recent papers 24,25 , we present a unified framework for the DM and 
MOND phenomena based on the rich and well-studied physics of superfluidity. The DM and 

“Consider a test particle orbiting a galaxy of mass Mb, in the low acceleration regime. Equating the centripetal 

acceleration v 2 /r to the MONDian acceleration ^/onOo = we obtain a velocity that is independent of 

distance, v 2 = \ZG^Mbao, in agreement with the flat rotation curves of spiral galaxies. Squaring this gives the 
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BTFR relation Mb = as an exact prediction. 














MOND components have a common origin, representing different phases of a single underlying 
substance. The central idea is that DM forms a superfluid inside galaxies, with a coherence 
length of galactic size. 

As is familiar from liquid helium, a superfluid at finite temperature (but below the critical 
temperature) is best described phenomenologically as a mixture of two fluids 26 > 27 > 28 : /) the 
superfluid, which by definition has vanishing viscosity and carries no entropy; ii) the “normal” 
component, comprised of massive particles, which is viscous and carries entropy. The fraction 
of particles in the condensate decreases with increasing temperature. Thus our framework 
naturally distinguishes between galaxies (where MOND is successful) and galaxy clusters (where 
MOND is not). Galaxy clusters have a higher velocity dispersion and correspondingly higher 
DM temperature. For rn ~ eV we will find that galaxies are almost entirely condensed, whereas 
galaxy clusters are either in a mixed phase or entirely in the normal phase. 

As a back-of-the-envelope calculation, we can estimate the condition for the onset of su¬ 
perfluidity ignoring interactions among DM particles. With this simplifying approximation, the 
requirement for superfluidity amounts to demanding that the de Broglie wavelength AdB ~ 1/ mv 
of DM particles should be larger than the interparticle separation t ~ (m/p) 1 / 3 . This implies 
an upper bound on the particle mass, m < (p/v 3 ) 1 / 4 . Substituting the value of v and p at 
virialization, given by standard collapse theory, this translates to 24,25 


m < 2.3 (1 + z vir ) 
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where M and z v j r are the mass and virialization redshift of the object. Hence light objects form 
a Bose-Einstein condensate (BEC) while heavy objects do not. 

Another requirement for Bose-Einstein condensation is that DM thermalize within galaxies. 
We assume that DM particles interact through contact repulsive interactions. Demanding that 
the interaction rate be larger than the galactic dynamical time places a lower bound on the 
interaction cross-section. For M = 10 12 /i _1 Mq and z v j r = 2, the result is 24,25 
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With rn < eV, this is just below the most recent constraint from galaxy cluster mergers 29 , 
though such constraints should be carefully reanalyzed in the superfluid context. 

Again ignoring interactions, the critical temperature for DM superfluidity is T c ~ mK, which 
intriguingly is comparable to known critical temperatures for cold atom gases, e.g., "Li atoms 
have T c ~ 0.2 mK. Cold atoms might provide more than just a useful analogy — in many 
ways, our DM component behaves exactly like cold atoms. In cold atom experiments, atoms are 
trapped using magnetic fields; in our case, DM particles are attracted in galaxies by gravity. 


3 Superfluid Phase 

Instead of behaving as individual collisionless particles, the DM is more aptly described as 
collective excitations: phonons and massive quasi-particles. Phonons, in particular, play a key 
role by mediating a long-range force between ordinary matter particles. As a result, a test 
particle orbiting the galaxy is subject to two forces: the (Newtonian) gravitational force and the 
phonon-mediated force. 

Specifically, it is well-known that the effective field theory (EFT) of superfluid phonon ex¬ 
citations at lowest order in derivatives is a P(X) theory 30 . Our postulate is that DM phonons 
are described by the non-relativistic MOND scalar action, 


(W ) 2 


P{X) ~ KXJ\X \; 


X = 6-m& 


2m 


(3) 



where A ~ meV to reproduce the MOND critical acceleration, and <J> is the gravitational po¬ 
tential. The fractional 3/2 power would be strange if Eq. (3) described a fundamental scalar 
held. As a theory of phonons, however, it is not uncommon to encounter fractional powers in 
cold atom systems. For instance, the Unitary Fermi Gas (UFG) 31,32 , which has generated much 
excitement recently in the cold atom community, describes a gas of cold fermionic atoms tuned 
at unitarity. The effective action for the UFG superfluid is uniquely fixed by 4d scale invariance 
at lowest-order in derivatives, £ufg(A 1) ~ X 5 / 2 , which is also non-analytic 33 . 

To mediate a force between ordinary matter, 9 must couple to the baryon density: 

£int = ~ a Yr ’ (^) 

Mpi 

where a is a dimensionless parameter. This term explicitly breaks the shift symmetry, but only 
at the 1/Mpi level and is therefore technically natural. From the superfluid perspective, Eq. (4) 
can arise if baryonic matter couple to the vortex sector of the superfluid, giving rise to operators 
~ cosdpb that preserve a discrete subgroup of the continuous shift symmetry 34,35,36 . 


3.1 Properties of the Condensate and Phonons 

The form of the phonon action uniquely fixes the properties of the condensate through standard 
thermodynamics arguments. At finite chemical potential, 6 = /it, ignoring phonon excitations 
and gravitational potential to zero, the pressure of the condensate is given as usual by the 
Lagrangian density, 

P(g) = |(2mg) 3 / 2 . (5) 

This is the grand canonical equation of state, P = P{p), for the condensate. Differentiating 
with respect to /i yields the number density of condensed particles: 

n = = A(2m) 3//2 /i 1//2 . (6) 

Combining these expressions and using the non-relativistic relation p = mn, we find 


P = 


12A 2 m 6 


(7) 


This is a polytropic equation of state P ~ p 1+1 / n with index n = 1/2. 

Including phonons excitations 9 = pt + 0, the quadratic action for is 


r A(2m.) 3 / 2 

£quad = 4 /x i / 2 




The sound speed can be immediately read off: 


c s = 
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3.2 Halo profile 


Assuming hydrostatic equilibrium, we can compute the density profile of a spherically-symmetric 
DM condensate halo: 
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Figure 2 Numerical solution of Lane-Emden equation, Eq. (11). 


Substituting the equation of state given by Eq. (7), and introducing the dimensionless variables 
p = po=i and r = y 32 nG^A J rrt s £> w ith po denoting the central density, Eq. (10) implies the 
Lane-Emden equation 

(£ 2 2')' = -« 2 = 1/2 . HD 

where ' = d/d£. The numerical solution, with boundary conditions 5(0) = 1 and 5/(0) = 0, is 
shown in Fig. 2. The superfluid density profile is cored, not surprisingly, and therefore avoids 
the cusp problem of CDM. 

The density is found to vanish at £i ~ 2.75, which defines the halo size: R = £i- 

Meanwhile the central density is related to the halo mass as 3 ' po = 4 ^ 3 > with 5/(£i) ~ 

—0.5. Combining these results, it is straightforward to solve for p q and R: 
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Remarkably, for m ~ eV and A ~ meV we obtain DM halos of realistic size! In the standard 
CDM picture a halo of mass Mdm = 10 12 M 0 has a virial radius of ~ 200 kpc. In our framework, 
the condensate radius can in principle be considerably smaller or larger depending on parameter 
values. For concreteness, in the remainder of the analysis we will choose as fiducial values 


m = 0.6 eV ; A = 0.2 meV . 

This implies a condensate radius of ~ 125 kpc for a halo of mass Mdm = 10 12 M 0 . 


(13) 


4 Phonon-Mediated MONDian Force 


Next we derive the phonon profile in galaxies, modeling the baryons as a static, spherically- 
symmetric localized source for simplicity. We first focus on the zero-temperature analysis, where 
the Lagrangian is given by the sum of Eqs. (3) and (4). In the static spherically-symmetric 
approximation, 9 = pt + <j>(r), the equation of motion reduces to 

V • (A^T V*) = , (14) 

where X(r) = p — m&(r) — - . This can be readily integrated: 
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There are two branches of solutions, depending on the sign of X. We focus on the MOND 
branch (with X < 0): 

<p' (r) = \fm ^JjL + \Jfi 2 + n 2 /m 2 ^j , 


where pi = p — m&. Indeed, for n/m 3> fi we have 

<t>'(r) ~ \Jn(r). 

In this limit the scalar acceleration on an ordinary matter particle is 
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To reproduce the MONDian result omond 
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we are therefore led to identify 


a 3 / 2 A = \JaoMp\ ~ 0.8 meV 
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which fixes a in terms of A through the critical acceleration. For the fiducial value A = 0.2 meV, 
we obtain a — 2.5. 

As it stands, however, the X < 0 solution is unstable. It leads to unphysical halos, with 
growing DM density profiles 24,25 . The instability can be seen by expanding Eq. (3) to quadratic 
order in phonon perturbations p = <p — (p(r), 
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The kinetic term ip 2 has the wrong sign for X < 0. (The X > 0 branch, meanwhile, is stable 
but does not admit a MOND regime 24,25 .) 

Since the DM condensate in actual galactic halos has non-zero temperature, however, we 
expect that the zero-temperature Lagrangian (Eq. (3)) to receive finite-temperature corrections 
in galaxies. At finite sub-critical temperature, the system is described phenomenologically by 
Landau’s two-fluid model: an admixture of a superfluid component and a normal component. 
The finite-temperature effective Lagrangian is a function of three scalars 38 : = E(A, B , Y). 

The scalar X, already defined in Eq. (3), describes the phonon excitations. The remaining scalars 
are defined in terms of the three Lagrangian coordinates ip 1 (x,t), I = 1, 2, 3 of the normal fluid: 


B = \J det ; 

Y = (d^6 + m5^j — m ~ p — m<3? + 4> + v ■ V0, (21) 

where u M eijRdaip 1 dpip J d^'ip K is the unit 4-velocity vector, and in the last step 

for Y we have taken the non-relativistic limit ~ (1 — $,17). By construction, these scalars 
respect the internal symmetries: i ) ip 1 —» ip 1 + c 1 (translations); ii ) ip 1 —>• R I j'ip J (rotations); 
in) ip 1 —>• ^(ip), with det = 1 (volume-preserving reparametrizations). 

There is much freedom in specifying finite-temperature operators that stabilize the MOND 
profile. The simplest possibility is to supplement Eq. (3) with the two-derivative operator 

AC = M 2 Y 2 = M 2 (fi + <P) 2 , (22) 


where we have specialized to the rest frame of the normal fluid, v = (1. This leaves the static 
profile given by Eq. (16) unchanged, but modifies the quadratic Lagrangian by M 2 p 2 , restoring 
stability for sufficiently large M. Specifically this is the case for 

1/4 


M > 


Am 3 / 2 


0.5 


10 11 M P 


o 


M h 


f a y/ 2 

\meV ) 


10 kpc 


1/2 


m. 


(23) 



which, remarkably, is of order eV! Hence, for quite natural values of M, this two-derivative 
operator can restore stability. Furthermore, this operator gives a contribution A P = M 2 /j, 2 
to the condensate pressure, which obliterates the unwanted growth in the DM density profile. 
Instead, the pressure is positive far from the baryons, resulting in localized, finite-mass halos 24,25 . 

5 Observational Implications 

We conclude with some astrophysical implications of our DM superfluid. 

Gravitational Lensing: In TeVeS 19 the complete absence of DM requires introducing a time-like 
vector field A M , as well as a complicated coupling between </>, A M and baryons in order to re¬ 
produce lensing observations. In our case, there is no need to introduce an extra vector, as the 
normal fluid already provides a time-like vector u^. Moreover, our DM contributes to lensing, 
so we are free to generalize the TeVeS coupling 24,25 . 

Vortices: When spun faster than a critical velocity, a superfluid develops vortices. The typical 
angular velocity of halos is well above critical 24,25 , giving rise to an array of DM vortices per¬ 
meating the disc 39 . It will be interesting to see whether these vortices can be detected through 
substructure lensing, e.g., with ALMA 40 . 

Galaxy mergers: A key difference with ACDM is the merger rate of galaxies. Applying Landau’s 
criterion, we find two possible outcomes. If the infall velocity Uj n f is less than the phonon sound 
speed c s (of order the viral velocity 24,25 ), then halos will pass through each other with negligible 
dissipation, resulting in multiple encounters and a longer merger time. If c; n f > c s , however, 
the encounter will excite DM particles out of the condensate, resulting in dynamical friction and 
rapid merger. 

Bullet Cluster: For merging galaxy clusters, the outcome also depends on the relative fraction 
of superfluid vs normal components in the clusters. For subsonic mergers, the superfluid cores 
should pass through each other with negligible friction (consistent with the Bullet Cluster), while 
the normal components should be slowed down by self interactions. Remarkably this picture 
is consistent with the lensing map of the Abell 520 “train wreck” 4 L42,43,44^ w hi c j 1 s how lensing 
peaks coincident with galaxies (superfluid components), as well as peaks coincident with the 
X-ray luminosity peaks (normal components). 

Dark-bright solitons: Galaxies in the process of merging should exhibit interference patterns (so- 
called dark-bright solitons) that have been observed in BECs counterflowing at super-critical 
velocities 45 . This can potentially offer an alternative mechanism to generate the spectacular 
shells seen around elliptical galaxies 46 . 

Globular clusters: Globular clusters are well-known to contain negligible amount of DM, and 
as such pose a problem for MOND 47 . In our case the presence of a significant DM component 
is necessary for MOND. If whatever mechanism responsible for DM removal in ACDM is also 
effective here, our model would predict DM-free (and hence MOND-free) globular clusters. 
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